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Abstract 

The mid-rapidity {dapN/dy at y=0) and total (upAr) production cross sections of 
J/f) mesons are measured in proton-nucleus interactions. Data collected by the 
HERA-B experiment in interactions of 920 GeV/c protons with carbon, titanium 
and tungsten targets are used for this analysis. The J/fi mesons are reconstructed 
by their decay into lepton pairs. The total production cross section obtained is 
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^'pN ~ ± 74 ± 46 nb/nucleon. In addition, our result is compared with previous 

measurements. 


PACS; 13.20Gd, 13.85Ni, 24.85+p 


1 Introduction 


Since the spectacular discovery of the J/ip particle [1,2] heavy-quarkonium 
production in hadron-hadron interactions has been in the focus of interest 
because it provides important information on both perturbative and non- 
perturbative aspects of Quantum Chromodynamics (QCD). Especially the 
unexpectedly large cross sections for J/ijj and t/’(2S') at large transverse mo¬ 
menta observed by the CDF experiment [3] renewed this interest and led to the 
development of the non-relativistic QCD (NRQCD) approach [4], which ex¬ 
tends the color-singlet model by including color-octet contributions [5,6,7]. A 
different approach to charmonium production is based on the color evaporation 
model [8,9]. Both models predict the energy dependence oi J/ip production. 

These theoretical predictions have to be compared to the experimental re¬ 
sults obtained in proton induced reactions. However, despite the large interest 
in this held, the experimental situation is far from being satisfactory. Cross 
section measurements at comparable energies differ well outside the quoted 
uncertainties. Thus a new measurement with low systematic errors is of inter¬ 
est. 

In this paper, we report on a measurement of the cross section for J/ijj produc¬ 
tion in interactions of 920 GeV/c protons with nuclei of atomic weight A using 
the HERA-B detector. The corresponding center of mass energy of the proton 
nucleon interaction is ^/s=4:1.6 GeV. A data sample of 210 million interactions 
on carbon, titanium and tungsten targets was recorded using a minimum bias 
trigger in the 2002-2003 HERA running period. The J/ijj mesons are detected 
in the inclusive reaction 

pA^ J/-ip X with J/'ll: e~^e~ or . 


The advantage of such a minimum bias data sample is the low systematic 
error due to the large (> 97%) trigger efficiency, the large angular coverage of 
HERA-B and its large reconstruction efficiency. On the other hand, the size 
of the sample is small, since charmonium cross sections are 4 — 5 orders of 
magnitude less than the total inelastic cross section. 
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2 Apparatus 


The HERA-B fixed-target spectrometer operated at the 920 GeV/c proton 
beam of the HERA storage ring at DESY using one or more wire targets in¬ 
serted into the beam halo. The detector was equipped with a vertex detector 
and extensive tracking and particle identification systems. It had a large geo¬ 
metrical coverage from 15 mrad to 220 mrad in the bending (horizontal) plane 
and 15 mrad to 160 mrad in the non-bending (vertical) plane. Fig. 1 shows a 
plan view of the detector in the configuration of the 2002-2003 data run. 

The target system [10] consisted of two stations of four wires each. The wires 
were positioned above, below, and on either side of the beam and were made 
from various materials including carbon, titanium and tungsten. The stations 
were separated by 40 mm along the beam direction. The wires were positioned 
individually in the halo of the stored proton beam and the interaction rate for 
each inserted wire was adjusted independently. 

The Vertex Detector System (VDS) [11] was a forward micro-strip vertex 
detector integrated into the HERA proton ring. It provided a precise mea¬ 
surement of primary and secondary vertices. The VDS consisted of 7 stations 
(with 4 stereo views each) of double-sided silicon strip detectors (50 x 70 mm, 
50 /im pitch) integrated into a Roman pot system inside a vacuum vessel and 
operated as near as 10 mm from the beam. An additional station was mounted 
immediately downstream of the 3 mm thick aluminum window of the vacuum 
vessel. 

The first station of the main tracker was placed upstream of the 2.13 Tm 
spectrometer dipole magnet. The remaining 6 tracking stations extended from 
the downstream end of the magnet to the electromagnetic calorimeter (ECAL) 
located 13 m downstream of the target. Each tracking station was divided 
into inner and outer detectors. The region starting from the beam pipe and 
extending up to 200 mm was covered by micro-strip gas chambers with GEM 
foils (inner tracker [ 12 ]) which, however, were not used in this analysis due 
to insufficient stability. The region outside the inner tracker was covered by 
the large area outer tracker (OTR) [13] consisting of ~95,000 channels of 
honeycomb drift cells. 

Particle identification was performed by a Ring Imaging Gherenkov detector 
(RIGH) [14], an electromagnetic calorimeter (EGAL) [15] and a muon de¬ 
tector (MUON) [16]. The RIGH used G 4 F 10 as radiator gas and two large 
spherical mirrors to project Gherenkov photons on the photon detector em¬ 
ploying multi-anode photomultipliers. The EGAL was based on “shashlik” 
sampling calorimeter technology, consisting of scintillator layers sandwiched 
between metal absorbers. It was subdivided in three different sections with 
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increasing cell size. In the radially innermost section, tnngsten was used as 
an absorber, and lead was used everywhere else. The MUON detector was 
segmented into four super-layers. Iron and concrete shielding extended from 
just behind the ECAL to the last MUON super-layer, except for gaps for the 
super-layers themselves. The hrst two super-layers consisted of three layers 
of tube chambers with different stereo angles. The last two super-layers each 
consisted of one layer of tube chambers with additional cathode pad readout. 



Figure 1. Plan view (bending plane) of the HERA-B detector. 


3 Data Sample and Trigger 


The present analysis is performed on single wire carbon, titanium and tungsten 
runs taken under stable conditions with a minimum bias trigger with a total 
of 182 million interactions (Table 1). The trigger required at least 20 hits in 
the RICH detector (compared to an average of 33 for a full ring from a /? = 1 
particle [14]) or an energy deposit of at least 1 GeV in the electromagnetic 
calorimeter and was sensitive to etrigger > 97% of the total inelastic cross 
section ainei- 


The integrated luminosity was determined [17] from the number of inelastic 
interactions Nimi using the expression C = Ninei/{(-trigger ■ o'inei) ■ The data were 
recorded at a moderate interaction rate of about 1.5 MHz which corresponds 
to 0.17 interactions per hlled bunch crossing. Therefore only about 10% of the 
events contain more than one interaction. The high data acquisition rate of 
about 1000 Hz allowed to record the bulk of the data within two weeks. 

The information delivered by the vertex detector, outer tracker, electromag¬ 
netic calorimeter and muon detector entered into this analysis. 
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Target 

Events 

C[yb-^) 

A(£)/£ 

C Below 1 

68.8 M 

277. 

4.8% 

C Inner 2 

20.5 M 

98.0 

4.8% 

Ti Below 2 

24.7 M 

30.9 

5.0% 

W Inner I 

67.6 M 

35.7 

3.8% 


Table 1 

Summary of minimum bias statistics. Target material, wire position with respect to 
the proton beam (inner wire corresponds to +X and outer wire to -X, see Fig. 1), 
integrated luminosities as well as the sample-specific luminosity errors are given. In 
addition, a 2.0% scaling error on luminosity has to be taken into account. 

4 Data Analysis 

4-1 Event Selection 


The search for J/ijj candidates is performed by analyzing the invariant mass 
spectrum of unlike-sign lepton pairs. The phase-space region of the J/'^, in 
rapidity y, covered by our measurement is —1.25 < y < -1-0.35. The corre¬ 
sponding range in Feynman xp is —0.225 < xp < -1-0.075. Throughout this 
paper, the rapidity is dehned in the proton-nucleon center-of-mass system. 
Soft cuts (see Table 2) are applied on the number of hits in the tracking sys¬ 
tem to select properly reconstructed tracks. Cutting on the dilepton vertex £t 
probability ensures that the two tracks have a common vertex. A cut on the 
transverse momentum of the leptons helps to reduce the background with a 
small loss (< 3%) in signal efficiency. 

The signal to background ratio can be improved by applying more stringent 
cuts on the lepton identification. In order to find the best set of cut values the 
signal significance S/\^S + B is maximized. S is the number of J/'^ surviving 
the selection cuts in the J/'^ Monte Carlo sample (see Sec. 5) scaled to obtain 
a comparable number of reconstructed J/'ip in Monte Carlo as in real data. 
The background B is evaluated from the like sign invariant mass spectrum in 
a 2.5 standard deviation window around the J/ip position. Such a procedure 
allows to adjust the lepton identification criteria in an unbiased way. 

Muon candidates are selected by requiring that the muon likelihood (nor¬ 
malized to one) derived from the muon hit information is greater than 0.2. 

In the case of electrons, the background is higher due to pions interacting 
in the ECAL and charged hadrons overlapping with energy depositions by 
neutral particles. We therefore require that the track position extrapolated 
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p^p 

e“’'e 

Number of hits in YDS 

> 6 

(12 hits on average) 



Number of hits in OTR 

> 10 

(40 hits on average) 



Transverse mom. pT (GeV/c) 

> 0.7 

Dilepton rapidity 

-1.25 <y < -hO.35 

Dilepton vertex fit probability 

> 10-5 

Muon likelihood (L^) 

> 0.2 


VDS-ECAL match 


|Ay| < 0.75 cell 
width 

E/p 


-2.0 < {E/p - 
0.98)/cj < -F3.5 


Table 2 

Event selection requirements. The lepton identification cuts (last three rows) are 
hxed by the optimization procedure described in the text. 

to the ECAL in the non-bending direction Y matches a cluster with |AY| 
< 0.75 cell width. Moreover, the ratio E/p, where E is the energy deposited 
in the electromagnetic calorimeter and p the momentum measured by the 
tracking system, provides good discrimination between electrons and hadrons. 
However, the electrons can emit Bremsstrahlung photons (BR) when passing 
the detector planes. BR photons emitted before the magnet are searched for 
in the ECAL and the momentum of the electron is corrected if such a cluster 
is found. Details of this procedure are explained in [18]. 

Additional energy losses result in a small deviation of the ratio E/p from 
unity. Consequently, a track is identified as an electron candidate if —2.0 < 
{E/p — 0.98)/(T < -1-3.5 with a ^ 6%. 


4-2 J/ip pP^p, 


The dimuon invariant mass distribution that results from combining all data¬ 
sets and applying the selection criteria discussed above is plotted in Fig. 2a. 
This distribution is fitted by a Gaussian plus a tail in the signal region and 
an exponential for the background. The tail takes into account the radiative 
process P/ip ^ [19] whereas the background shape is suggested by the 

mass spectrum of like-sign lepton pairs (see Fig. 2b). 






Figure 2. a) Unlike sign pairs fitted with a Gaussian plus a radiative tail for the 
signal and an exponential for the background; b) like sign dimuon candidates fitted 
with an exponential. 

A fit with this function results in 100 ±12 J/'if) ^ iji~ decays. To indicate 
the influence of the radiative tail, the mass distribution is htted with a simple 
Gaussian resulting in 94 ± 12 events. 

The peak position and the width (FWHM) are 3.093 ± 0.005 GeV/c^ and 
85 ± 10 MeV/c^, respectively. The width is in good agreement with the width 
expected from the J/'ip Monte Garlo simulation (see Sec. 5). These £t results 
are used as input to the £t of the individual signals obtained from the three 
different data samples shown in Fig. 3a-c. The sum of the three signals (see 
Table 3) is in good agreement with the £t to the total sample. 


J/lp 

total sample 

subsample 

Fit including radiative tail 

100±12 


carbon 


25±6 

titanium 


16±4 

tungsten 


58±9 

J/'ip —^ e+e~ 

total sample 

subsample 

Fit including radiative tail 

57±13 


carbon 


32±8 

titanium 


2±4 

tungsten 


24±10 


Table 3 

Signal yield for J/ip and J/ip —> e'''e 
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Figure 3. Invariant mass of unlike sign muon pairs from a) carbon, b) titanium, and 
c) tungsten sample. 

4.3 J/ip e+e“ 


In contrast to muons, the momentum measurement of electrons is distorted by 
energy losses due to the emission of Bremsstrahlung in the detector planes be¬ 
fore the electromagnetic calorimeter. Consequently, the dielectron mass distri¬ 
bution shows a more pronounced tail towards smaller masses. The description 
of this tail taking into account Bremsstrahlung losses as well as the radia¬ 
tive tail due to the decay J/'0 —»• e'''e “7 was adjusted by using our dilepton- 
triggered J/'ifj ^ e'''e“ sample [20] which was recorded in the same data taking 
period with the same detector set-up. Fig. 4a displays the dielectron invariant 
mass distribution after applying the selection criteria listed in Table 2. Using 
an exponential to describe the background we obtain 57 ± 13 J/'^ —> e'''e“ 
events (Table 3). As for the muons, the exponential background can be moti¬ 
vated by the like sign mass spectrum shown in Fig. 4b. 
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Figure 4. a) Invariant mass of unlike sign electron pairs fitted with a Gaussian plus 
a radiative tail for the signal and an exponential for the background; b) like sign 
dielectron candidates fitted with an exponential. 

5 Efficiency Determination 


A Monte Carlo simulation is used to determine the J/'0 —*■ efficiencies. 

The Monte Carlo samples for pA J/ip + X are generated in two steps. 
First, a cc pair is generated with Pythia 5.7 [21] such that a J/'ip is always 
produced. The generated events are reweighted such that the resulting cross 
sections conform to the parameterizations 


da 

dpx'^ 


1 + 


f 35 ■ 71 ■ pt\ 
V256 -(pt); 


-6 


and 


da 

dxp 


{I-MY. 


The average transverse momentum of (px) = 1.29 ±0.01 GeV/c is taken from 
a preliminary HERA-B analysis of the J/p) sample recorded with a dilepton 
trigger [20]. The slight increase of the average transverse momentum with the 
mass of the nucleus is neglected. We take c = 6.38 ± 0.24 as measured in 
proton-silicon collisions at 800GeV/c [22]. The J/ip^s are generated without 
polarization. The influence of a possible J/ip polarization within limits given 
by the experimental results [22-25] is taken into account as a systematic error. 
After the generation of the J/'ip the remaining energy is given as input to 
Fritiof 7.02 [26] which generates the underlying event taking into account 
further interactions inside the nucleus. 
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The detector response is simulated with the Geant 3.21 package [27]. Realistic 
detector efficiencies, readout noise and dead channels are taken into account. 
The simulated events are processed by the same reconstruction codes as the 
data. The total efficiencies, including track reconstruction as well as the influ¬ 
ence of the selection cuts, are summarized in Table 4. The trigger efficiency 
for events containing & J/ip meson is above 99%. 


Target 


e(e+e ) 

C Below 1 

0.337 ± 0.003 

0.256 ± 0.003 

C Inner 2 

0.339 ± 0.004 

0.261 ± 0.004 

Ti Below 2 

0.325 ± 0.003 

0.247 ± 0.002 

W Inner 1 

0.317 ± 0.002 

0.235 ± 0.002 


Table 4 

Total efficiency of reconstruction and selection cuts in the rapidity interval —1.25 < 
y < -1-0.35. The uncertainties quoted are the statistical uncertainties of the MC 
simulation. 


6 Results 


6 .1 Combined cross sections per nucleus 


The first step in determining the production cross sections is to combine, for 
each target material, the results of the two different final states. This is possible 
because the two measurements are statistically independent and compatible. 
The weighted average takes into account the statistical errors as well as those 
systematic errors which depend on the lepton species (contributions 1 - 2 of 
Sec. 6.2). The visible cross sections per nucleus, i.e. the cross sections measured 
in the rapidity interval of —1.25 < y < -1-0.35, are given by 

^ _ 

Br{J/p, ^ £+r) ■ Ci ■ A 


where 

• Ni, Ci and Ci are the measured J/ip yield, the efficiency (Table 4) and the 
integrated luminosity (Table 1) for a particular target, 

• Br is the branching fraction which is (5.88 ± 0.10) % for J/ip ^ pCp,~ and 
(5.93 ± 0.10) % for J/pj e+e" [28]. 

The combined cross sections ^ summarized in Table 5, are almost inde¬ 

pendent of the assumptions made on the differential distributions (see Sec. 5). 
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However, for the total cross sections an extrapolation factor 1// must be ap¬ 
plied which depends on the shape of the rapidity distribution. This factor 
/ = 0.631 ± 0.010 is determined from the rapidity distribution of the J/ijj 
mesons generated in the Monte Carlo simulation. 

The total cross sections obtained by extrapolation of the visible cross sections 

^j/y_ 

^pA — 7 


are also given in Table 5. 


Target 

A 

l\apA (pb) 

apA (pb) 

carbon 

12.011 

4.0 ±0.7 ±0.3 

6.3 ±1.1 ±0.5 

titanium 

47.867 

17. ±6. ±1. 

27. ±9. ±2. 

tungsten 

183.84 

75. ±11. ±5. 

118. ± 18. ±8. 


J/'0 cross 
statistical 


Table 5 

Visible AapA (in the rapidity interval —1.25 < y < -1-0.35) and total apA 
sections per nucleus of atomic weight A. The errors quoted indicate the 
and systematic uncertainties, respectively. 


6.2 Systematic Uncertainties 


The total systematic uncertainty of the production cross sections is composed 
of the following contributions: 

(1) Particle identihcation cuts are varied over a wide range and, for each set 
of cut values, the cross section is evaluated. The variation of the cross 
section corresponds to 2.5% for J/ip jj,'^jj,~ and 2.0% for J/ip —> e’''e“, 
respectively. 

(2) Allowing for extreme variations of both signal shape and background 
parameterization, we estimate the rms uncertainty on signal counting by 
dividing the difference of the extreme values by vT2 resulting in 2.5% for 
J/pj ^ jjUij,~ and 8% for J/pj ^ e’''e“. 

(3) The luminosity uncertainties of each data sample are shown in Table 1. 

(4) The overall scaling error of the luminosity is 2.0% [17]. 

(5) The track reconstruction efficiency is tested separately for VDS and OTR 
using decays. One of the charged pions is reconstructed based on 
VDS (or OTR) and RICH/ECAL information only. Applying the same 
procedure to Monte Carlo data, we obtain an uncertainty of the track 
reconstruction efficiency of 1.5% per track. 


13 




(6) The uncertainty of the Monte Carlo production model within our kine- 
matical range is 2.5%. This includes variations of the px shape as well as 
of the polarization parameter. Varying the exponent c of the xp distribu¬ 
tion within errors has negligible inflnence on onr reconstrnction efficiency. 
Redncing the exponent c from 6.38 to 5.0 corresponds to a 1% change in 
the efficiency. 

(7) Each of the branching fractions for J/'ip and J/p) —> e’''e“ has 

an nncertainty of 1.7% [28]. 

(8) The nncertainty on the extrapolation to the fnll rapidity range is deter¬ 
mined by varying the exponent of the xp parameterization within the 
errors to be 1.5%. 

(9) Using as parameterization of the A-dependence, the 

nncertainty of o; = 0.96 ± 0.01 [29] corresponds to a 3.6% nncertainty of 
the cross section per nncleon. 


6.3 Cross sections per nucleon 


As a last step, the resnlts of the three data samples are combined to determine 
the cross section per nncleon The cross section for J/p) prodnction on a 

nnclear target of atomic weight A is parameterized as 


n, 


j/y 

pA 


= a. 


Jli) 

pN 




By htting this expression to onr data we can determine the cross section per 
nucleon. However, the limited statistics of our sample does not allow for a pre¬ 
cise measurement of the A-dependence. Therefore we £x the parameter to the 
mid-rapidity measnrement a = 0.96 ± 0.01 of E866 [29]. Fig. 5 demonstrates 
that onr resnlt is in good agreement with this assnmption [x^/NDF = 1.5/2). 
For this fit only the statistical and the relevant systematic nncertainties (1 - 
3 of Sec. 6.2) are taken into acconnt. 

As a resnlt we obtain 

a'pj^ = 663 ± 74 ± 46 nb/nncleon 


for the total J/p) prodnction cross section per nncleon. The systematic nncer¬ 
tainty qnoted comprises all contribntions described in Sec. 6.2. 

The mid-rapidity cross section can be derived by applying a factor p = 1.50 ± 
0.03 which corrects for the difference between the differential cross section 
at y = 0 and the visible cross section, which is an average of the interval 
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J/\|i total cross section 



Figure 5. The total J/i!) cross section on nuclear targets C, Ti and W. The error 
bars indicate the statistical and systematic uncertainties 1-3 explained in Sec. 6.2. 
The line shows the result of the fit described in the text using a = 0.96. 

—1.25 < y < +0.35 with Ay = 1.6 


da. 


j/ip 

pN 


dy 


y=0 


Aa^''^ 

Ay 


■ V 


392 ±44 + 27 nb/nucleon. 


7 Comparison with other results 


Comparing our results to the measurements of previous experiments in proton- 
induced interactions requires correcting all results for the same branching frac¬ 
tions and A-dependence. All results shown in Fig. 6 are updated for the latest 
values of the branching fractions quoted in Sec. 6 and a target mass depen¬ 
dence assuming a = 0.96 + 0.01. The systematic uncertainties are recalculated 
accordingly. 

It is obvious that the experimental results are far from being consistent. Cross 
sections measured at nearby beam energies are not compatible when taking 
into account just the errors quoted by the experiments. Also shown on the 
figure is a fit to the published data on proton-induced J/tp and 'ip{2S) produc¬ 
tion in the context of a next-to-leading order NRQCD calculation [50]. Several 
fits are performed with different subsets of measurements and different parton 
distribution functions (PDF). The best description of the data is obtained 
with the PDF MRST2002 [51]. 

At a/s =41.6 GeV, the fit gives a cross section of 

a^l^ = 502 ± 44 nb/nucleon 
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J/y cross section da/dy|y^o 


* pp 


C 

o 

O 10'^r O pA 
c = • HERA-B 

!q I 


10 " 


■D 

■D 


10 



10 


10 


Vs (GeV) 


J/y total cross section 



Figure 6. J/'0 production cross sections in proton-induced interactions [22-24, 
30-47]. pp and pA measurements are indicated by different symbols, a) Differen¬ 
tial cross section dapN/dy at y=0] b) total cross section cPpAr together with the fit 
described in the text. 


which is in reasonable agreement with the result presented here. The error 
takes into account both fit uncertainties as well as systematics due to the 
parton distribution function. Details of the procedure together with tables of 
all experimental input used are the subject of a separate paper [52]. 
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8 Summary 


A data sample taken with a minimum bias trigger by the HERA-B experiment 
in interactions of 920 GeV/c protons with carbon, titanium and tungsten tar¬ 
gets has been used to determine the production cross section oi J/ip mesons. 
In our data sample we hud 100 ±12 J/ip ^ and 57 ± 13 J/ip ^ e+e“ 
candidates. After correcting for the efficiency of the selection criteria within 
the range in rapidity of —1.25 < y < ±0.35 and combining both hnal states, 
we obtain the visible cross sections per nucleus within the detector acceptance. 
Using ■ ylo.96±o.oi parameterize the A-dependence, and extrapo¬ 

lating this measurement to the full y range, the total production cross section 
per nucleon is 

— 663 ± 74 ± 46 nb/nucleon. 


For the cross section at mid-rapidity we obtain 


da. 


pN 


dy 


j/=0 


392 ± 44 ± 27 nb/nucleon. 


An NRQCD based £t of published results on J/ip and ip{2S) production pre¬ 
dicts = 502 ± 44 nb/nucleon at 920 GeV/c, in reasonable agreement 

with this measurement. 
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